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Abstract 



CP asymmetry in B — > 4>Ks decay is studied in the special context of su- 
persymmetry theories, in which the charginos play an important role. We 
find that in addition to the gluino, the chargino can also make large contribu- 
tions to CP asymmetry in B —* <j)Ks decay. After considering the constraints 
from B — > J/ipKs decay, we study three special scenarios: (a). Large mixing 
between left-handed charm and top squarks(LL mixing); (b). Large mix- 
ing between right-handed charm and top squarks (RR mixing); (c). Large 
mixing between both left-handed charm-top squark and right-handed charm- 
top squark(LL+RR mixing). We show quantitatively that because of large 
squark mixing within the second and third generations, an 0(1) effect on CP 
violation in B — > <j)Ks is possible. 

PACS numbers: 13.25.Hw; 12.60.Jv; 11.30.Pb 
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I. INTRODUCTION 



In standard model, all the CP asymmetries in the B system can be accounted for with 
a single phase in Cabbibo-Kobayashi-Maskawa mixing matrix [1]. It predicts that the time 
dependent CP asymmetry in B — > 4>K$ decay is the same as the CP asymmetry in B — > 
J/vpKs decay. That is: 

Sj/^Ks = S^s = sin2p = 0.731 ± 0.056(51/) 

Any difference between these two decays would be evidence of new physics beyond the 
standard model. Finding CP violation in the B system is a very important goal in B 
factories. 

BaBar [2] and Belle [3] measurements on CP asymmetries in B — > J/ipK s decay have 
provided the first evidence for CP violation in the B system. 

sin2(3(B -> J/ipK s ) = 0.741 ± 0.067 ± 0.033 (BaBar) 

sin2f3(B -> J/^Ks) = 0.733 ± 0.057 ± 0.028 (Belle) ' ' 1 

The average of their results is: 

S J/i!Ks = sin2(3 = 0.734 ± 0.054, 

This result is in good agreement with the SM prediction [4]. Therefore, we can conclude 
that CP is significantly violated and the CKM phase angle is the dominant source of CP 
violation. 

However, the CKM mechanism is challenged by the recent measurements of the BaBar 
and Belle collaborations in B — > (j)K s decay. 

sin2(3(B -> <pK s ) = +0.45 ± 0.43 ± 0.07 (BaBar) 

sin2P(B -> 4>K S ) = -0.96 ± 0.50lg;?? (Belle) ^ ' 

The average of these two measurements implies: 

S^s = sin2(3 = -0.15 ± 0.33. 

which deviates from the SM prediction by 2.7c. The 2a difference between Sj/^k s an d S^k s 
is expected to be an indication of new physics. Many papers have explored potential new 
physics contributions [5-9] to B — > <pK s decay. Among them, supersymmetry is one of the 
most significant candidates to account for this 2a deviation. 

In B — > J/ipKs decay, standard model contributions dominate the tree level diagrams 
and any new physics enters only in loop diagrams, which are suppressed naturally. As a 
result, in order to have a significant supersymmetric contribution, a large SUSY CP violat- 
ing phase is required. Such large phases usually do not exist in most of the supersymmetric 
models. Therefore the CP asymmetry in Sj/^k s is dominated by the standard model contri- 
bution. Unlike the case of B — > J/ipKs, there are no tree level diagrams in B — > <pKs decay 
because b — > sss , the dominant process in B — > <f)Ks, is induced only at one loop level 
in the standard model. Both SM and SUSY contributions appear at one loop level. Thus 
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SUSY contributions are more significant in B — > 4>Ks decay and are expected to explain the 
2c difference. 

SUSY contributions to B — > decay come mainly from penguin and box diagrams 
containing gluino, chargino and neutralino loops. The contributions from neutralino loops 
are much smaller than the contributions from gluino and chargino loops. We ignore contribu- 
tions from neutralinos in our computation. There are several papers investigating on SUSY 
contributions to B — > <pKs decay [10-18], especially the gluino loop contributions [19-24]. 
Only a few include the chargino diagrams, which can also affect the result substantially [25]. 
The main purpose of this paper is to study chargino contributions to the CP asymmetry in 
B — > 4>Ks decay. The chargino contributions depend on the details of the model. However, 
we try to be as general as possible in our analysis. To emphasize the chargino loop contribu- 
tions, we work in a special basis where the down-type squark mass matrix is diagonalized. 
There is no intergeneration mixing in the down-type squark sector. In this scenario, gluino 
contributions to B — > <pKs decay are ruled out and chargino contributions are the dominant 
contributions. 

This paper is organized as follows: In the next section, we present the formula for CP 
violation in the B system and introduce the effective field theory for the b — > s transition. 
Section III is devoted to the discussion of supersymmetric contributions to the CP asym- 
metry in B — > (j)K s process. Section IV includes the conclusion and the loop functions are 
showed in Appendix. 



II. CP VIOLATION INB-> <f)K s DECAY 

In the B system, the time dependent CP asymmetry in B — > <pKs decay is defined as 

T(B°(t) - <f)Ks) - T(B°(t) - <j)Ks) 



l <t>K S 



r(B°(t) - <pK s ) + r(5°(t) - (^Ks) 



Csk, cos AM B .t + Ssk, sin AM B .t 



(3) 



where AMs d is the mass difference between the two neutral B mesons. C<$,k s and S$k s 



are 



_ |A(0^ 5 )| 2 -1 _2/m 



5 K<i>Ks 



|A(0K 5 )| 2 + 1' 
p and q are the mixing parameters defined as: 

\B ± >=p\B° > ±q\B° > 
\(j,K s is a ratio of decay amplitudes: 



|A(0^ 5 )| 2 + 1 



(4) 



A 



4,Kg 



A(B" -> <PK S ) 
A(B° -> <f)Ks) 



(5) 



From above formulae, we notice that there are two sources of CP violation [26]. One is 
from the mixing ||| ^ 1, the other is from the ratio of decay amplitudes |^| ^ 1 [4,27]. In 
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the SM, loop diagrams with W boson and top quark dominate the b — > sss decay channel. 
In this channel, A^ s = a 6m {b°^4,Ks) = y*vf = 1- There is no CP violation from the ratio 
of decay amplitudes. Thus, CP asymmetry in B — > (j)Kg decay comes mainly from mixing , 
which is the same with B — > J/ipK s decay. 

Supersymmetry can contribute to mixing, AB = 2 transitions, and to decay AB = 1 
transitions. SUSY contributions to the decay amplitudes are through penguin and box 
diagrams [28], while contributions to the mixing parameters are only through box diagrams. 
To concentrate on the chargino contributions, we choose a basis, in which the down-type 
squark mass is diagonalized. There is no intergeneration mixing for down-type sector. The 
gluino contributions are ruled out by the choice of basis. SUSY contributions are mainly 
chargino contributions. 

The effective Hamiltonian in AB = 1 transitions is: 

H = j:[C^)0^)}+C g 9 (6) 
i=i 

As in Ref. [23,29], the operators we choose are: 

1 = {s alll P L c p ){c^fP L b a \ 
C 2 = {s alt ,P L c a ){c^P L b p ), 

3 = {s al »P L b a )(s pl »P L sp), 

04 = (s a ^P L bp)(s^P L S a ), 
06 = {SalfiPLb^ispYPRSa), 

O g = £- 2 m b (s a T^ u P R b,)G^ (7) 

where a and (3 are color indices, Pl — (1 — 7s)/2, and a^ v = |[7 /1 , 7"]- The contribution 
from the operator 1 = ^m^Sad ^ u PRb a ) F^ u is very small and is ignored in this paper. 
Operators proportional to s-quark mass m s are also ignored because of the small s-quark 
mass. 

In the SM, the Wilson coefficients 6*3^6 m P> 4>Ks decay have been computed to NLL 
order, and C g has been computed to LL order. The values we used are in ref. [30]. 

In SUSY, up-squark masses are taken to be non-degenerate. We use the mass eigenbasis 
in the computation. There are large mixing between the second and third generations in 
up-type sector. The mass matrix takes the form: 



/ 
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m Lct 




















m lt 








m t A t cot(3 















m Ru 
























m Rc 


m Rct 




V 








mtAtcot/3 





m*L 


m Rt 


J 



(8) 



This mass matrix can be diagonalized by a unitary matrix Ru. That is: 
Rim\Ru = diag(ml 1 ,ml 2 ,ml 3 ,m% 1 ,m% 2 ,m 2 m ), 
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It is useful to define the 6x3 matrices such that 



(r^ I ')aj \Ru)aii (J^^Jai \-R"u)ai+3- 

where a = 1, ... ,6 denotes the mass eigenstates and i = u,c,t labels the gauge eigenstates. 
The chargino mass matrix is 

M 2 V2M w sinf3 \ 
V2M w cosp n )' ^ ' 

where we assume the mass of the W-ino field, M 2 , is real and positive. The non observation 
of the neutron electric dipole moment (EDM) strongly constrains the phase of \i. A non- zero 
phase of /i does not significantly affect our results, we keep the phase of /i to be zero for 
sake of simplicity. 

The chargino mass matrix can be diagonalized by a biunitary transformation 

diag(m~±, m-±) = UM^±V T 

where U and V are unitary matrices defined as 

jj _ ( cosOjj sindjj 
\ —sinOu cosOu 





y = ( cos ^v sinOy \ 
I —sinOy cosOy J 

with mixing angles [29] 

2M W [M 2 +/j, 2 ± (Mf - n 2 ) cos 2(3 + 2/iM 2 sin 2(3} 1 ' 2 



sin 29 uv 



m 2 ± — m 2 ± 

Xi X2 



M 2 - f i 2 T 2M 2 v cos2[3 

COs26'[/ i y = 2 \ . (12) 

mz± — mi± 

Xi X2 

We calculated Wilson coefficients of SUSY diagrams in Fig. 1 with chargino loops 
[23,31-33]. 

2 

/^iSUSY \ " ®"w l s~i*iks u*iks\ ( /~iihs Tjihs\ I /^i*jhs Tj*jhs\(/^jkb Tjjkb\ 

6 3 =2^ -J^T^ G u L ~ H u R ){ G u L ~ H u R A^u L ~ H u R )\y J u L - H t R ) 

i,h,j,k *" l Xj 

C(a; m - fe t m-. , ^m fife ,m & . i x rriu k ,m Sj ) 

+ E - K 3 R k s ){Gt b - K j R kb )Ci{x m ^ 3 ) 

j,k Xj 



\ ^ 2ol w (X\ , „*jk s TT*jks\/ (-t*jkb TT*jkb\ 

+ Qm 2 \ U L A% n u R ) 

j,k ym Xj 
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c 4 = l^-^r( G u L ~ H u J R )i G u L ~ H u R )^i{x m . k 

j,k Xj 

a 2 

/^iSUSY \ "> "ui / /~i*iks Tj*iks\ Tjihs Tj*jhs ( /~<jkb Tjjkb\ 

C 5 =2^ Ui ) H u L H u L V^u L - H t R ) 

i,h,j,k Xj 

®-s®-w / f~t*jks TJ*jks\( /-i*jkb Tj*ikb\/~i I 

j,k Xj 

_|_ 2,(X w (Xi ^jxjfcg jj*jks^^y*jkb Jj*jkb^ 

j,k 9m Xj 



/^iSUSY \ ^ Oi s (X w , ^y*jks Tj*jks\f^i*jkb zj*jkb\sy / \ 

°6 ~ L^-7r-?r\^u L - n u R ){^u L ~ n u R )^i{x m - k , m - ) 

j,k Zm Xj 

/^<SUSY \ ^ OZ w 7l , ^^jks TT*iks\( /-t*ikb zj*jkb\ r-> / \ 

°9 - I^O^T^ul ~ U u R ){^u L ~ U u R )^l{X m - k ,m s ) 

j,k MU Xj 

+ E^(C - H*J R ks WJ R kb -f D 2 {x V] ) (13) 

j,k Zl7l Xj m b 

(14) 



where 

m^, j — 1,2 are chargino masses, and x mA)ms = m 2 A /m 2 B . The loop functions are given 
in Appendix. 

In addition to Wilson coefficients, the hadronic matrix elements of the operators are also 
needed in computing the SUSY contributions to the CP asymmetry in B — > <f>K s decay. We 
use the results in ref. [23,34] 

(<PK s \O 3A \B d ) = \H (l + ^ 
{ct>Ks\O b \B" d ) = hi 
(<PK s \O 6 \B d ) = 

(<j>K S \o 9 m = ^ e ^lw- (16) 

(17) 

where the analytical expression for H is in Appendix B of ref. [23]. The value of k is 
approximately -1.1 [20,23]. 
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III. NUMERICAL RESULTS 



With the formulae in section II, we can numerically analyze the chargino contributions 
in B system. There are three kinds of mixing in the up-squark mass matrix: Left-Left, 
Right-Right and Left-Right. Each mixing contributes a mixing angle and a phase angle in 
Scf>K s computation. However, the phase of /i is strongly constrained by the neutron electric 
dipole moment (EDM). We have checked that our results have no substantially dependence 
on fi. To simplify our computation, we take the phase of /i to be zero in this paper. 

The agreement between the CP asymmetry in B — > J/ipKs decay observed by BaBar 
and Belle with the SM prediction leaves very little room for new physics contributions in 
Sj/^k s - It strongly constrains the SUSY contributions to Sj/^k s and strongly suppresses 
the CP violating phase in B — > 4>K S system [35], which will be noticed in the following 
analysis. 

In order to simplify the analysis, we consider three special cases. First, left-left mixing 
dominates, LL scenario. Right-right mixing is very small and can be ignored in this scenario. 
The second scenario is right-right mixing dominant, RR scenario, while the last one has 
comparable left-left and right-right mixing contributions, LL + RR scenario. We need to 
consider both contributions in the last case. Because of the large parameter space, we only 
show some results for RR + LL scenario with specific simplifications. 



A. Left-Left Mixing Dominant 

In this scenario, left-left mixing is the only sizable mixing and there is only one phase 
and one mixing angle. We take the mixing matrix to be: 

( cos6 23 sin0 23 e i<t>L \ ( . 

\ -sin6 23 cos0 23 e-^ )' { } 

where 823 is the mixing angle. The subscript 23 represents the left-left mixing and is the 
mixing between second and third rows in up-squark mass matrix. 0l is the phase angle. We 
expect that these mixing and phase angles can significantly alter the value of S^k- We have 
checked that the precise values of M 2 and /i do not have a substantial effect on our results. 
We use a range of values of M 2 and /x from 2 GeV to about 400 GeV. We obtain almost the 
same results. Thus, we take M 2 and /1 to be 200 GeV in most of our calculations. 
In this scenario, we have parameters 

m L 2 » m U ' m X2 > m Xl » tan P > #23 , 0L • 

with mi 3 < mi 2 and < m^ 2 . As pointed in ref . [35] , the imaginary parts of this matrix 
are stringently constrained by B — > J/v/jK s decay. In order for the mixing angle 6 23 to be less 
constrained, the phase 4>l we use in our computation is very severe constrained, (pi < 10 -2 . 

There are seven parameters. Their effects on S^Ks are different. In Fig. 2, we plot S^k 
as a function of mi and within the range —0.7 < S^k < 0.7. We fixed tanf3 = 5, 
and mi 2 = 5 TeV. 9 23 is taken to be 7r/4, and (pL = 0.01. In frame (a), The value of k is 
equal to -1.1 and m^ 2 = lOm^. In frame (b), k has the same value as in frame (a), but 
m^ 2 = lOOm^. We change k to -2 in frame (c) and frame (d) with m^ 2 = 10m^ 2 in frame 
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(c) and rriy, 2 = 100m^ 2 in frame (d). Within the shaded region, S^k is larger than -0.7, but 
smaller than 0.7. From this graph, we notice that: 

• There are many assumptions in the process of k estimation, k = —1.1 is just an 
approximate value. We present results for different values of k to show the dependence 
of S^Ks 011 K - The main SUSY contributions in this scenario come from penguin 
diagrams. The hadronic matrix of the corresponding operator is proportional to k and 
k surely has some effect on S^Ks- This can be observed through the difference among 
graphs in Fig. 2. 

• We show how large the effect of the parameter k is. We take k to be -2 in frame (c) 
and frame (d). Except for these two frames, k is fixed to -1.1 in all other graphs. 

• We focus on chargino loop contributions. Any changes in the chargino parameters 
have a large influence on our results. As is seen in frame (c) and frame (d), increasing 
m^ 2 changes the allowed values significantly. 

• There are big gaps around the line = m^. We refer the reader to the loop 
functions in Eq. 14, which are presented explicitly in Appendix. When the chargino 
and squark masses are equal, the functions in the Appendix take a different form and 
cause the big gaps in the graphs. 

• The small discontinuity in frame (a) and frame (b) at = 500 GeV arises for the 
same reason. That is because m^ 2 is ten times bigger than m^, and at the point of 

= 500 GeV, m^ 2 = 5 TeV, which is equal to m^. 

We further explore the dependence of S^k s on the phase angle 4>l and the mixing angle 
6>23 in Fig. 3 and Fig. 4. We reduce the value of 4>l to 0.005 in Fig. 3. We continue to take 
tan/? = 5, mi 2 = 5 TeV and m^ 2 = lOm^. As shown in the plot, with a smaller <p L angle 
the allowed region is suppressed < 200 GeV. In Fig. 4, # 2 3 is reduced to tt/6 and 4>l 
is fixed to be 0.01. The allowed region changes as a result. Fig. 4 presents an even bigger 
change as 9 decreases from 7r/4 to tt/6. 

There is one parameter left: tan/3, which is very important. From Eq. 12 and Eq. 15, we 
can see that tan/3 affects not only the mixing angles, but also the Wilson Coefficients. Any 
change in tan/3 alters S^Ks significantly. We plot Fig. 5 to explore the dependence of S^Ks 
on tan/3. Apart from tan/3, all the other parameters remain the same as the parameters in 
frame (a) of Fig. 2. In frame (a), we take tan/3 to be 20, while in frame (b), tan/3 is set to 
be 40. Several features are worth noting: 

• By comparing two graphs in this figure, we find that the allowed region permits larger 
masses with the increase in tan/3. 

• As for graph in frame (b), we just plot S^k s i n a region with the masses of m~ Lz and 
m^j smaller than 1.5 TeV. There is a large allowed region beyond 1.5 TeV, which does 
not appear in this graph. Still, we can see the strong dependence of S$k s on tan/3. 
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B. Right-Right Mixing Dominant 



This corresponds to the case where left-left mixing is very small and can be ignored. The 
dominant mixing is right-right mixing, the mixing matrix is: 

( cos6 56 sin$ 56 e^ R \ ( , 

\ -sin6 m cos9 m e- i(t> K J ' 1 ' 

The free parameter are now 

m R 2 i m R 3 i m X2, m Xi > tan P> #56 , <t>R- 

with the same constraints m^ g < and < m^ 2 . 

Unlike the Left-Left dominant scenario, this time, we fix m^ 3 = 300 GeV. In Fig. 6, 
we plot S^Ks as a function of m^, 2 and m^. tan/3 is still 5. We try different values of 
(f) within the range <fi < 0.01, and find that in almost the whole range, the value of S^Ks 
is approximately 0.73, which is the same as the SM prediction. SUSY contributions are 
very small. We relax the constraint and take to be tt/4. 6> 56 is also 7r/4. In frame 
(a), rriy. 2 = lOm^, while in frame (b), m^ 2 is 20 times bigger than m^. From this graph, 
we notice that, in contrast to the left-left mixing case, the shaded region in frame (b) is 
concentrated toward small masses side than that in frame (a) . As shown in Eq. 12, the 
difference between the chargino masses, m^ 2 — affects the mixing angles in the chargino 
sector. The mixing angles determine the values of S^Ks- We allow m^ 2 to be proportional 
to m^j such that the larger the mass of , the bigger the mass difference of m^ 2 and . 
It is the fact that m^ 2 in frame (b) is ten times bigger than the corresponding m^ 2 in frame 

(a) make the allowed S^Ks region in frame (a) concentrate on larger , while in frame 

(b) , is relatively small. 

At = 300 GeV, we observe a small peak in the graphs. As in the left-left mixing 
case, at this point equals to , which uses a difference loop function. 

We also wish to explore the dependence of S^,k s on other parameters. Unlike in left-left 
mixing case, this time instead of studying the dependence of S^k s on each parameter, we 
show a graph scanning the whole parameter space: 

m L 2 ' m U ' m xi » m X2 » #23, <Pr- 
Within the constraints m^ 2 > and m^ 2 > m^, we take: 

300 GeV < m k2 < 2500 GeV (20) 
250 GeV < m k3 < 1000 GeV 
100 GeV < rrifr < 500 GeV 
150 GeV < m^ 2 < 1000 GeV 

< # 23 < tt/2 

< 4> R < 2tt 

The results are present in Fig. 7 with S^Ks vs m x,2- 111 frame (a), the tan/5 is 15, and in 
frame (b) it is increased to 50. In frame (a), we notice that many of the points are around 
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0.73, the value of standard model prediction. This behavior tells us with the parameters in 
Eq. 20 SUSY contributions are usually small and cannot alter S^Ks substantially. What we 
are interested in is the regions of parameter space where SUSY contributions are dominant. 
In these regions, the value of S^Ks can be as low as negative -0.98. As tan/5 increases, SUSY 
contributions are increased. As shown in frame (b), when tan/3 is increased to 50, more 
points on the scatter plot deviate from the SM prediction. The mixing angles 9jj and 9y 
in the chargino sector vary with tan/3 and increase the SUSY contributions to S^Ks- There 
are very few points within m^ 2 < 300 GeV. That is because we require m^ 2 > m^, which 
limits the probability of small m^ 2 . 

C. Left-Left Mixing + Right-Right Mixing 

Left-left mixing and right-right mixing are comparable in this scenario. They both 
make contributions to S^k s ■ There are two mixing angles, two phase angles and four mass 
parameters. To simplify the computation, we assume that left and right squarks have the 
same masses in same generation, that is: m^ 2 = m^ 2 = rfi2 and m^ 3 = mjj 3 = m 3 . With 
these assumptions, the mixing angles 9l in the left-left mixing and 9r in the right-right 
mixing are the same, 9l = 9 R = 9 and also the phase angle, 0l = 0_r = 0. Therefore, the 
free parameters in this scenario are reduced to: 

rh 2 , m 3 , m^ 2 , m^, tan/3, 9, 0. 

with the same constraints rh 3 < m 2 and < m^ 2 . 

We display our main results in a few plots. In Fig. 8, we fixed rfi2 = 5 TeV and the mixing 
angle 9 is set to be 7r/4. The phase angle is 0.01. In frame (a), we use m^ 2 = 10m^ 1 
, while in frame (b) m^ 2 is 50 times bigger than m^. This graph is very similar to the 
corresponding graph in left-left mixing scenario, Fig 2, except here we have a larger allowed 
parameter space. As in the LL dominant scenario, chargino masses affect S$k s significantly. 
This is shown in the difference between graphs in frame (a) and frame (b). As expected, left- 
left mixing provides the dominant contribution and right-right mixing adds to the allowed 
region. The big gaps around line rh 3 = and the small gap at rh 2 = m^ 2 are much more 
obvious than in left-left dominant scenario. We can conclude that left-left plus right-right 
mixing together make a big contribution to S^k s - 

S^Ks a l so depends on the mixing and phase angles. Fig. 9 shows the effect of phase 
angle 0, while Fig. 10 shows the effect of the mixing angle 9. We keep rfi2 = 5 TeV and 
m X2 = IOttt,^. In Fig. 9, the phase angle is reduced to 0.005 and mixing angle is fixed 
to be 7r/4, while in Fig. 10, 9 changes to n/6 and is kept at 0.01. As in the LL case, as 
decreases the allowed parameter space decreases and as 9 becomes smaller the allowed 
masses also become smaller. 

Fig. 11 is the scanned graph. We scan the same parameter space as in the right-right 
mixing case. We find that in frame (a), although there is still a concentration of points 
around the standard model value, more points deviate from their SM prediction than with 
only right-right mixing. This effect is even more obvious in frame (b). This also shows a 
stronger contribution in this scenario. 



10 



IV. CONCLUSION 



We have studied the supersymmetric contributions, especially the chargino loop contri- 
butions, to CP asymmetry in B — > <f)Ks decay. To emphasize chargino loop contributions, 
we work in a special basis where the down-type squark mass matrix is diagonal. In this 
special basis, gluino loop contributions, which are very important in most SUSY frames, are 
ruled out. Chargino loops are the main contributions to S^Ks- 

Within the allowed parameter region, we studied three special cases. Left-left mixing 
dominant, right-right mixing dominant and left-left plus right-right mixing dominant. In 
the LL case, chargino contributions to the CP asymmetry can be much larger than the 
SM contributions in some regions of parameter space. The CP asymmetry S^k s varies 
significantly as the chargino masses vary. In the RR case, things are very similar with the 
LL case, except the allowed region is totally different from the LL case. The LL plus RR 
mixing gives the largest contribution to S^k s - 

We consider SUSY contributions to CP asymmetry S^Ks with only chargino loops. In 
most SUSY scenarios, gluino loops will provide very large contributions, which may be larger 
than chargino contributions for some parameters. Gluino, chargino and even neutralino loops 
make up the whole SUSY contribution. 

b — > S7 decay puts very strong limits on the allowed parameter space. It excludes some 
regions in our graphs. The calculation of b — * 37 is very model dependent. We are not 
including the b — * S7 constraints in our paper. But we can still conclude that SUSY is a 
candidate to explain the CP asymmetry deviation in B — > (pKs and B — > J/ipKs and the 
contributions of chargino loops play an important role in some regions of SUSY parameter 
space. 
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APPENDIX A: LOOP FUNCTIONS 




A 2 log A B 2 \ogB 



(A - 1)(A - B)(A - C) (B - 1)(B - A)(B - C) 
C 2 \ogC 



(Al) 



(C-1)(C- A)(C-B)' 
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A log A 



B\ogB 



(A - l)(A - B)(A - C) (B - 1)(B - A)(B - C) 

ClogC 

(C- 1)(C- A){C-BY 



2x 3 - 9x 2 + 18x - 11 - 6 log a; 

36(1 - a;) 4 ' 

-lQx 3 + 45x 2 - 36a; + 7 + Qx 2 (2x - 3) log a; 
36(1 - a;) 4 

—x 3 + 6a; 2 — 3a; — 2 — 6a; log x 
6(1 -a;) 4 ' 

—x 2 + 1 + 2x log x 
{x- l) 3 ' 
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FIG. 1. Box and penguin diagrams in b — ► sss. 
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FIG. 2. In left-left mixing scenario, S^Ks i n the plane of rfiL3 y s ^xi- ^ ne shaded region 
corresponds to —0.7 < S^k s < 0-7 with <fi = 0.01 and = tt/4. (a), k = —1.1, m X 2 = 10m x i; (b). 
k = —1.1, m X 2 = 100m x i; (c). k = —2., m X 2 = 10m x i; (d). k = —2., m X 2 = 100m x i; 
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FIG. 3. In left-left mixing scenario, S^k s m the plane of rfiL3 vs rh Xl . The shaded region 
corresponds to —0.7 < S^k s < 0.7 with <fi = 0.005, 6 = ir/4 and k = —1.1. m X 2 = 10m x i 
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FIG. 4. In left-left mixing scenario, S^,k s 111 the plane of fh,L3 vs fh Xl . The shaded region 
corresponds to —0.7 < S^k s < 0-7 with <p = 0.01, 8 = it/6 and k = —1.1. m X 2 = 10m x i 
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= 20; (b) tan/3 = 40. 
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FIG. 8. In LL + RR mixing case. S^,k s hi the plane of ttils vs rh xi . The shaded region 
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FIG. 9. LL + RR case. S^Ks in the plane of ttilz vs m Xl . The shaded region corresponds to 
-0.7 < S^Ks < 0.7 with <f> = 0.005, 9 = vr/4 and k = -1.1. fn x2 = 10m x i 
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FIG. 10. LL + RR case. S^Ks in the plane of friL3 vs rh Xl . The shaded region corresponds to 
-0.7 < S^Ks < 0-7 with ^ = 0.01, 6 = tt/6 and k = — 1.1. m X 2 = 10m x i 
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FIG. 11. LL + RR case. The CP asymmetry S^k s scanned within parameter space over the 
region mentioned in the text. We show the results in the plane of S^Ks vs ^ls with (a), tan/3 = 
15 and (b). tan/3 = 50 
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